A weaponized phage suppresses competitors in historical and modern metapopulations of pathogenic bacteria.

Bacteriophages, the viruses of bacteria, are proposed to drive bacterial population dynamics, yet direct evidence of their impact on natural populations is limited. Here we identified viral sequences in a metapopulation of wild plant-associated Pseudomonas spp. genomes. We discovered that the most abundant viral cluster does not encode an intact phage but instead encodes a tailocin - a phage-derived element that bacteria use to kill competitors for interbacterial warfare. Each pathogenic Pseudomonas sp. strain carries one of a few distinct tailocin variants, which target variable polysaccharides in the outer membrane of co-occurring pathogenic strains. Analysis of historic herbarium samples from the last 170 years revealed that the same tailocin and receptor variants have persisted in the Pseudomonas populations for at least two centuries, suggesting the continued use of a defined set of tailocin haplotypes and receptors. These results indicate that tailocin genetic diversity can be mined to develop targeted "tailocin cocktails" for microbial control. One-Sentence Summary: Bacterial pathogens in a host-associated metapopulation use a repurposed prophage to kill their competitors.

Richness and density jointly determine context dependence in bacterial interactions.

Pairwise interactions are often used to predict features of complex microbial communities due to the challenge of measuring multi-species interactions in high dimensional contexts. This assumes that interactions are unaffected by community context. Here, we used synthetic bacterial communities to investigate that assumption by observing how interactions varied across contexts. Interactions were most often weakly negative and showed a phylogenetic signal among genera. Community richness and total density emerged as strong predictors of interaction strength and contributed to an attenuation of interactions as richness increased. Population level and per-capita measures of interactions both displayed such attenuation, suggesting factors beyond systematic changes in population size were involved; namely, changes to the interactions themselves. Nevertheless, pairwise interactions retained some explanatory power across contexts, provided those contexts were not substantially divergent in richness. These results suggest that understanding the emergent properties of microbial interactions can improve our ability to predict the features of microbial communities.

Microbe-associated molecular pattern recognition receptors have little effect on endophytic Arabidopsis thaliana microbiome assembly in the field.

Plant microbiome structure affects plant health and productivity. A limited subset of environmental microbes successfully establishes within plant tissues, but the forces underlying this selectivity remain poorly characterized. Transmembrane pattern recognition receptors (PRRs), used by plants to detect microbe-associated molecular patterns (MAMPs), are strong candidates for achieving this selectivity because PRRs can potentially interact with many members of the microbiome. Indeed, MAMPs found in many microbial taxa, including beneficials and commensals, can instigate a robust immune response that affects microbial growth. Surprisingly, we found that MAMP-detecting PRRs have little effect on endophytic bacterial and fungal microbiome structure in the field. We compared the microbiomes of four PRR knockout lines of Arabidopsis thaliana to wild-type plants in multiple tissue types over several developmental stages and detected only subtle shifts in fungal, but not bacterial, ß-diversity in one of the four PRR mutants. In one developmental stage, lore mutants had slightly altered fungal ß-diversity, indicating that LORE may be involved in plant-fungal interactions in addition to its known role in detecting certain bacterial lipids. No other effects of PRRs on α-diversity, microbiome variability, within-individual homogeneity, or microbial load were found. The general lack of effect suggests that individual MAMP-detecting PRRs are not critical in shaping the endophytic plant microbiome. Rather, we suggest that MAMP-detecting PRRs must either act in concert and/or are individually maintained through pleiotropic effects or interactions with coevolved mutualists or pathogens. Although unexpected, these results offer insights into the role of MAMP-detecting PRRs in plant-microbe interactions and help direct future efforts to uncover host genetic elements that control plant microbiome assembly.

Iron drives anabolic metabolism through active histone demethylation and mTORC1.

All eukaryotic cells require a minimal iron threshold to sustain anabolic metabolism. However, the mechanisms by which cells sense iron to regulate anabolic processes are unclear. Here we report a previously undescribed eukaryotic pathway for iron sensing in which molecular iron is required to sustain active histone demethylation and maintain the expression of critical components of the pro-anabolic mTORC1 pathway. Specifically, we identify the iron-binding histone-demethylase KDM3B as an intrinsic iron sensor that regulates mTORC1 activity by demethylating H3K9me2 at enhancers of a high-affinity leucine transporter, LAT3, and RPTOR. By directly suppressing leucine availability and RAPTOR levels, iron deficiency supersedes other nutrient inputs into mTORC1. This process occurs in vivo and is not an indirect effect by canonical iron-utilizing pathways. Because ancestral eukaryotes share homologues of KDMs and mTORC1 core components, this pathway probably pre-dated the emergence of the other kingdom-specific nutrient sensors for mTORC1.

Dynamic genetic adaptation of Bacteroides thetaiotaomicron during murine gut colonization.

To understand how a bacterium ultimately succeeds or fails in adapting to a new host, it is essential to assess the temporal dynamics of its fitness over the course of colonization. Here, we introduce a human-derived commensal organism, Bacteroides thetaiotaomicron (Bt), into the guts of germ-free mice to determine whether and how the genetic requirements for colonization shift over time. Combining a high-throughput functional genetics assay and transcriptomics, we find that gene usage changes drastically during the first days of colonization, shifting from high expression of amino acid biosynthesis genes to broad upregulation of diverse polysaccharide utilization loci. Within the first week, metabolism becomes centered around utilization of a predominant dietary oligosaccharide, and these changes are largely sustained through 6 weeks of colonization. Spontaneous mutations in wild-type Bt also evolve around this locus. These findings highlight the importance of considering temporal colonization dynamics in developing more effective microbiome-based therapies.

Evolutionary implications of host genetic control for engineering beneficial microbiomes.

Engineering new functions in the microbiome requires understanding how host genetic control and microbe-microbe interactions shape the microbiome. One key genetic mechanism underlying host control is the immune system. The immune system can promote stability in the composition of the microbiome by reshaping the ecological dynamics of its members, but the degree of stability will depend on the interplay between ecological context, immune system development, and higher-order microbe-microbe interactions. The eco-evolutionary interplay affecting composition and stability should inform the strategies used to engineer new functions in the microbiome. We conclude with recent methodological developments that provide an important path forward for both engineering new functionality in the microbiome and broadly understanding how ecological interactions shape evolutionary processes in complex biological systems.

The design of mapping populations: Impacts of geographic scale on genetic architecture and mapping efficacy for defense and immunity.

Genome-wide association studies (GWAS) have yielded tremendous insight into the genetic architecture of trait variation. However, the collections of loci they uncover are far from exhaustive. As many of the complicating factors that confound or limit the efficacy of GWAS are exaggerated over broad geographic scales, a shift toward more analyses using mapping panels sampled from narrow geographic localities ("local" populations) could provide novel, complementary insights. Here, we present an overview of the major complicating factors, review mounting evidence from genomic analyses that these factors are pervasive, and synthesize theoretical and empirical evidence for the power of GWAS in local populations.

Plant genetic effects on microbial hubs impact host fitness in repeated field trials.

Although complex interactions between hosts and microbial associates are increasingly well documented, we still know little about how and why hosts shape microbial communities in nature. In addition, host genetic effects on microbial communities vary widely depending on the environment, obscuring conclusions about which microbes are impacted and which plant functions are important. We characterized the leaf microbiota of 200 Arabidopsis thaliana genotypes in eight field experiments and detected consistent host effects on specific, broadly distributed microbial species (operational taxonomic unit [OTUs]). Host genetic effects disproportionately influenced central ecological hubs, with heritability of particular OTUs declining with their distance from the nearest hub within the microbial network. These host effects could reflect either OTUs preferentially associating with specific genotypes or differential microbial success within them. Host genetics associated with microbial hubs explained over 10% of the variation in lifetime seed production among host genotypes across sites and years. We successfully cultured one of these microbial hubs and demonstrated its growth-promoting effects on plants in sterile conditions. Finally, genome-wide association mapping identified many putatively causal genes with small effects on the relative abundance of microbial hubs across sites and years, and these genes were enriched for those involved in the synthesis of specialized metabolites, auxins, and the immune system. Using untargeted metabolomics, we corroborate the consistent association between variation in specialized metabolites and microbial hubs across field sites. Together, our results reveal that host genetic variation impacts the microbial communities in consistent ways across environments and that these effects contribute to fitness variation among host genotypes.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity.

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection.

Gene duplication is increasingly recognized as an important mechanism for the origination of new genes, as revealed by comparative genomic analysis. However, how new duplicate genes contribute to phenotypic evolution remains largely unknown, especially in plants. Here, we identified the new gene EXOV, derived from a partial gene duplication of its parental gene EXOVL in Arabidopsis thaliana. EXOV is a species-specific gene that originated within the last 3.5 million years and shows strong signals of positive selection. Unexpectedly, RNA-sequencing analyses revealed that, despite its young age, EXOV has acquired many novel direct and indirect interactions in which the parental gene does not engage. This observation is consistent with the high, selection-driven substitution rate of its encoded protein, in contrast to the slowly evolving EXOVL, suggesting an important role for EXOV in phenotypic evolution. We observed significant differentiation of morphological changes for all phenotypes assessed in genome-edited and T-DNA insertional single mutants and in double T-DNA insertion mutants in EXOV and EXOVL. We discovered a substantial divergence of phenotypic effects by principal component analyses, suggesting neofunctionalization of the new gene. These results reveal a young gene that plays critical roles in biological processes that underlie morphological evolution in A. thaliana.

Functional biology in its natural context: A search for emergent simplicity.

The immeasurable complexity at every level of biological organization creates a daunting task for understanding biological function. Here, we highlight the risks of stripping it away at the outset and discuss a possible path toward arriving at emergent simplicity of understanding while still embracing the ever-changing complexity of biotic interactions that we see in nature.

Assessing the potential to harness the microbiome through plant genetics.

Microbial communities are influenced by a complex system of host effects, including traits involved in physical barriers, immunity, hormones, metabolisms and nutrient homeostasis. Variation of host control within species is governed by many genes of small effect and is sensitive to biotic and abiotic environments. On the flip side, these host impacts seem targeted on particular microbial species, with that impact percolating through the microbial community. There is not yet evidence that the nature and strength of these interactions differs between fungal and bacterial communities, or among different compartments of the plant. The challenge of deciphering how systems of host traits impact systems of microbial associates is vast but holds promise for developing novel strategies to improve plant health.

Metabolic Profile Discriminates and Predicts Arabidopsis Susceptibility to Virus under Field Conditions.

As obligatory parasites, plant viruses alter host cellular metabolism. There is a lack of information on the variability of virus-induced metabolic responses among genetically diverse plants in a natural context with daily changing conditions. To decipher the metabolic landscape of plant-virus interactions in a natural setting, twenty-six and ten accessions of Arabidopsis thaliana were inoculated with Turnip mosaic virus (TuMV), in two field experiments over 2 years. The accessions were measured for viral accumulation, above-ground biomass, targeted and untargeted metabolic profiles. The phenotypes of the accessions ranged from susceptibility to resistance. Susceptible and resistant accessions were shown to have different metabolic routes after inoculation. Susceptible genotypes accumulate primary and secondary metabolites upon infection, at the cost of hindered growth. Twenty-one metabolic signatures significantly accumulated in resistant accessions whereas they maintained their growth as mock-inoculated plants without biomass penalty. Metabolic content was demonstrated to discriminate and be highly predictive of the susceptibility of inoculated Arabidopsis. This study is the first to describe the metabolic landscape of plant-virus interactions in a natural setting and its predictive link to susceptibility. It provides new insights on plant-virus interactions. In this undomesticated species and in ecologically realistic conditions, growth and resistance are in a permanent conversation.

Natural Bacterial Assemblages in Arabidopsis thaliana Tissues Become More Distinguishable and Diverse during Host Development.

To study the spatial and temporal dynamics of bacterial colonization under field conditions, we planted and sampled Arabidopsis thaliana during 2 years at two Michigan sites and surveyed colonists by sequencing 16S rRNA gene amplicons. Mosaic and dynamic assemblages revealed the plant as a patchwork of tissue habitats that differentiated with age. Although assemblages primarily varied between roots and shoots, amplicon sequence variants (ASVs) also differentiated phyllosphere tissues. Increasing assemblage diversity indicated that variants dispersed more widely over time, decreasing the importance of stochastic variation in early colonization relative to tissue differences. As tissues underwent developmental transitions, the root and phyllosphere assemblages became more distinct. This pattern was driven by common variants rather than those restricted to a particular tissue or transiently present at one developmental stage. Patterns also depended critically on fine phylogenetic resolution: when ASVs were grouped at coarse taxonomic levels, their associations with host tissue and age weakened. Thus, the observed spatial and temporal variation in colonization depended upon bacterial traits that were not broadly shared at the family level. Some colonists were consistently more successful at entering specific tissues, as evidenced by their repeatable spatial prevalence distributions across sites and years. However, these variants did not overtake plant assemblages, which instead became more even over time. Together, these results suggested that the increasing effect of tissue type was related to colonization bottlenecks for specific ASVs rather than to their ability to dominate other colonists once established.IMPORTANCE Developing synthetic microbial communities that can increase plant yield or deter pathogens requires basic research on several fronts, including the efficiency with which microbes colonize plant tissues, how plant genes shape the microbiome, and the microbe-microbe interactions involved in community assembly. Findings on each of these fronts depend upon the spatial and temporal scales at which plant microbiomes are surveyed. In our study, phyllosphere tissues housed increasingly distinct microbial assemblages as plants aged, indicating that plants can be considered collections of tissue habitats in which microbial colonists-natural or synthetic-are established with differing success. Relationships between host genes and community diversity might vary depending on when samples are collected, given that assemblages grew more diverse as plants aged. Both spatial and temporal trends weakened when colonists were grouped by family, suggesting that functional rather than taxonomic profiling will be necessary to understand the basis for differences in colonization success.

The study of host-microbiome (co)evolution across levels of selection.

Microorganismal diversity can be explained in large part by selection imposed from both the abiotic and biotic environments, including-in the case of host-associated microbiomes-interactions with eukaryotes. As such, the diversity of host-associated microbiomes can be usefully studied across a variety of scales: within a single host over time, among host genotypes within a population, between populations and among host species. A plethora of recent studies across these scales and across diverse systems are: (i) exemplifying the importance of the host genetics in shaping microbiome composition; (ii) uncovering the role of the microbiome in shaping key host phenotypes; and (iii) highlighting the dynamic nature of the microbiome. They have also raised a critical question: do these complex associations fit within our existing understanding of evolution and coevolution, or do these often intimate and seemingly cross-generational interactions follow novel evolutionary rules from those previously identified? Herein, we describe the known importance of (co)evolution in host-microbiome systems, placing the existing data within extant frameworks that have been developed over decades of study, and ask whether there are unique properties of host-microbiome systems that require a paradigm shift. By examining when and how selection can act on the host and its microbiome as a unit (termed, the holobiont), we find that the existing conceptual framework, which focuses on individuals, as well as interactions among individuals and groups, is generally well suited for understanding (co)evolutionary change in these intimate assemblages. This article is part of the theme issue 'The role of the microbiome in host evolution'.

Population Genetics of the Highly Polymorphic RPP8 Gene Family.

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

Characterizing both bacteria and fungi improves understanding of the Arabidopsis root microbiome.

Roots provide plants mineral nutrients and stability in soil; while doing so, they come into contact with diverse soil microbes that affect plant health and productivity. Despite their ecological and agricultural relevance, the factors that shape the root microbiome remain poorly understood. We grew a worldwide panel of replicated Arabidopsis thaliana accessions outdoors and over winter to characterize their root-microbial communities. Although studies of the root microbiome tend to focus on bacteria, we found evidence that fungi have a strong influence on the structure of the root microbiome. Moreover, host effects appear to have a stronger influence on plant-fungal communities than plant-bacterial communities. Mapping the host genes that affect microbiome traits identified a priori candidate genes with roles in plant immunity; the root microbiome also appears to be strongly affected by genes that impact root and root hair development. Our results suggest that future analyses of the root microbiome should focus on multiple kingdoms, and that the root microbiome is shaped not only by genes involved in defense, but also by genes involved in plant form and physiology.

Genome-wide association studies on the phyllosphere microbiome: Embracing complexity in host-microbe interactions.

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe-related traits were developed in the relatively simple context of binary host-microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome-wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant-microbe interactions will influence the choice of traits to be mapped.

Genome-wide association study reveals new loci involved in Arabidopsis thaliana and Turnip mosaic virus (TuMV) interactions in the field.

The genetic architecture of plant response to viruses has often been studied in model nonnatural pathosystems under controlled conditions. There is an urgent need to elucidate the genetic architecture of the response to viruses in a natural setting. A field experiment was performed in each of two years. In total, 317 Arabidopsis thaliana accessions were inoculated with its natural Turnip mosaic virus (TuMV). The accessions were phenotyped for viral accumulation, frequency of infected plants, stem length and symptoms. Genome-wide association mapping was performed. Arabidopsis thaliana exhibits extensive natural variation in its response to TuMV in the field. The underlying genetic architecture reveals a more quantitative picture than in controlled conditions. Ten genomic regions were consistently identified across the two years. RTM3 (Restricted TEV Movement 3) is a major candidate for the response to TuMV in the field. New candidate genes include Dead box helicase 1, a Tim Barrel domain protein and the eukaryotic translation initiation factor eIF3b. To our knowledge, this study is the first to report the genetic architecture of quantitative response of A. thaliana to a naturally occurring virus in a field environment, thereby highlighting relevant candidate genes involved in plant virus interactions in nature.

The molecular genetic basis of herbivory between butterflies and their host plants.

Interactions between herbivorous insects and their host plants are a central component of terrestrial food webs and a critical topic in agriculture, where a substantial fraction of potential crop yield is lost annually to pests. Important insights into plant-insect interactions have come from research on specific plant defences and insect detoxification mechanisms. Yet, much remains unknown about the molecular mechanisms that mediate plant-insect interactions. Here we use multiple genome-wide approaches to map the molecular basis of herbivory from both plant and insect perspectives, focusing on butterflies and their larval host plants. Parallel genome-wide association studies in the cabbage white butterfly, Pieris rapae, and its host plant, Arabidopsis thaliana, pinpointed a small number of butterfly and plant genes that influenced herbivory. These genes, along with much of the genome, were regulated in a dynamic way over the time course of the feeding interaction. Comparative analyses, including diverse butterfly/plant systems, showed a variety of genome-wide responses to herbivory, as well as a core set of highly conserved genes in butterflies as well as their host plants. These results greatly expand our understanding of the genomic causes and evolutionary consequences of ecological interactions across two of nature's most diverse taxa, butterflies and flowering plants.